Berry phase for a spin 1/2 in a classical fluctuating fleld 
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The effect of fluctuations in the classical control parameters on the Berry phase of a spin 1/2 
interacting with a adiabatically cyclically varying magnetic fleld is analyzed. It is explicitly shown 
that in the adiabatic limit dephasing is due to fluctuations of the dynamical phase. 
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Berry phase IT| and related geometrical phases 0, ll| 
have received renewed interest in recent years due to sev- 
eral proposal for their use in the implementation of q uan- 
tum computing gates 0, H H S H M US El El El El • 
Such interest is motivated by the belief that geometric 
quantum gates should exhibit an intrinsic fault tolerance 
in the presence of external noise. Such belief is based on 
the heuristic argument that being Berry phases geometri- 
cal in their nature, i.e. proportional to the area spanned 
in parameter space, any fluctuating perturbation of zero 
average should indeed average out. Although this argu- 
ment seems convincing to the best of our knowledge it 
has not been quantitatively probed so far. In particu- 
lar although several papers have investigated 
aspects of Berry phases in the presence of quantum ex- 
ternal noise we are not aware of any in which the effect 
of classical noise in a simple model of qubit, namely a 
spin 1/2 interacting with an external classical field with 
a fluctuating component has been analyzed. This is pre- 
cisely the aim of this paper. For such system the effects 
of classical fluctuations in the control parameter on both 
geometric and dynamic phases is studied and their im- 
pact on dephasing analyzed. Our system consists of a 
spin 1/2 in the presence of an external static magnetic 
field, whose Hamiltonian, in appropriate units, takes the 
form 

H{t) = \^{t)-a (1) 

where a — (ctj;, dj,, cr^), di are the Pauli operators 
and B(t) = Bo{t)h{t) with the unit vector n = 
(sini? cos sini9 sin (/?, COS!?). The classical field B(i) 
acts as an external control parameter, as its direction 
and magnitude can be experimentally changed. When 
varied adiabatically the instantaneous energy eigenstates 
follow the direction of h and therefore can be expressed 
as 

T„) - e-^'^/2cos^|t)+e^'^/'sin^U) 

lU = e-W2sin^|t)_e^W2eos^|i) (2) 

where | t), | i) are the eigenstates of the CTz operator. 

When the time evolution is cyclic i.e. when after a time 
T we have B(T) = B(0) the energy eigenstates acquire 



a phase factor which contains a geometric correction to 
the dynamic phase: 

ITn(T))=e''e^^-|T„(o)) (3) 

where the dynamic phase S = B(){t)dt and the Berry 
phase can be expressed in term of the so called Berry 
connection = |Va| Tn) as jb = / " where 
A is the set of control parameters. In our specific example 
A = An analogous expression holds for A^, rela- 

tive to the I in) state. It is straightforward to calculate 
the components of A: 

Al = -A'l = z(t„ \d/dip\ t„) - i cosz? (4) 

Al = -4 =j(t„ \d/dd\ U) =0 (5) 

It is important to note that while the eigenenergies de- 
pend on Bo{t) the eigenstates depend only on n{t). As 
a consequence the Berry phase depends only on z9, ip. A 
standard example is a slow precession of B at an angle 
around the z axis with angular velocity uj ~ 2t:/T <^ Bq. 
A straightforward calculation shows that 

7T = -7i = / A^df = TTCOsi? (6) 
Jo 

Note that the the Berry phase is proportional to the 
solid angle subtended by B with respect to the degener- 
acy B = 

We are now in the position to extend our analysis to 
the case in which the magnetic field contains a fluctuating 
component. In this case Hamiltonian JQ) is modified as 
follows 

Hit) = ^BT-a^^{B{t)+K{t)).a (7) 

where we have divided the total magnetic field By into 
an average component B experimentally under our con- 
trol and a fluctuating field K. We will analyze the case 
in which B is a field of constant amplitude which under- 
goes a cyclic evolution while the components of K are 
random processes with zero average and small amplitude 



2 



compared to B, in order to consider lowest order cor- 
rections. Finally we will assume that the fluctuations 
are characterized by timescales such that the adiabatic 
approximation holds. We will show that this is not an 
unphysical restriction. 




FIG. 1: Generic path of the total magnetic field Bt. The 
field B is under experimental control and undergoes a cyclic 
evolution while K is a small random fiuctuating field. Note 
that in general K fiuctuates both in the direction parallel and 
perpendicular to B. 

Corrections to © have a twofold origin as the fluc- 
tuating field K modifies both the connection A and the 
path. To first order correction the connection is 



di9 



■5d = 



= 2 (cosi?o 



(5i?sini?o) 



(8) 
(9) 



where iSq is the polar angle of the average field B while 
SiD is the fiuctuation of the polar angle due to the fiuc- 
tuating field K. In order to analyze the corrections to 
the path, with no loss of generality, we again consider 
the case of a slow precession of B around the z axis. In 
the presence of K the line element dX will have also a 
component perpendicular to the orbit of B . However as 
the connection A has zero component in the direction 
we can restrict our attention to the f component of dX. 
To this end we write 



dX = ipdt ^ [ifQ + 5ip)dt 



(10) 



In 

•^0 



l|10|) ipQ is the average angular velocity. In our case 
= 2t:/T while Sip is the first order correction due to 
K: in particular when K fluctuates in the same direction 
of B the precession speed increases while it decreases in 
the case of fluctuation in the opposite direction. Note 
that fluctuations in the path contain only corrections in 
Lp while fluctuations in the connection depend only on 
This is independent from any approximation but is due 
to the structure of the connection. 



We can now express the Berry phase in the presence 
of noise as 



IB — I 
10 



{A^i^o) + SA^){(po + Sip)dt 



(11) 
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where the average Berry phase 7^ coincides with 7^ in 
the absence of noise and it has been assumed 5(^(0) = 0. 
The last term in Hll|l is a non-cyclic contribution which 
appears when, due to the presence of K, Bt does not 
return to its original. In this case, instead of the geomet- 
rical phase deflnition given by Berry, which assumes that 
the Hamiltonian is periodic, we have to use the deflnition 
by Samuel and Bhandari |l8l | about non cyclic evolution. 
If this is done the third term does not appear and eq. Hll|) 
becomes: 



75=75-77;/ sindoSddt 
J Jo 



(12) 



In order to proceed a physical model for the noise is 
needed, in other words a stochastic process for K must 
be assigned. Given the probability distribution for the 
field it is straightforward to calculate the distribution for 
the Berry phase. 

As a first step we express the trigonometric functions 
appearing in (|12|l in terms of the fluctuating fleld com- 
ponents Ki. This will be useful in calculating the prob- 
ability distribution for 7^. Let j Vo bo polar 
angles of B; i9 , those of Bx; and 5d , Sip the flrst or- 
der differences between the polar angles of the two flelds. 
Moreover let S = |B| be the modulus of B. If we expand 
in Taylor series cos we obtain: 



cos(z?o + S-d) = cos-do - (5i5sin??o = ^ + ^ - • K 

B B B 



and therefore 



-5d sin -iJo 



rB K 



B 

Substituting equation H13|l in H12|l wc find: 



7B = 7b + 



T 



T r 



B 



S3 
S3 



B K 



dt 



(13) 



(14) 



From this expression it is possible to find the probability 
distribution for 73, once that for Ki is known. 

We will assume that the fluctuating field Ki is a 
Ornstein-Uhlenbeck process, i.e it is gaussian, stationary 
and markovian with a lorentzian spectrum whose band- 
width Ti we assume to be much less than the Bohr fre- 
quency of our energy eigenstates. However while we must 
have Lu <^ B and Ti ^ B we have no restriction on the 
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relative value of vs Fj. In order to allow for the possi- 
bility of anisotropic noise we assume that -K3 has variance 
(T3 and bandwidth T3 while Ki and K2 have ai2 and ri2. 



With these assumption we found that the distribution 
for 7b is a gaussian whose mean value, as said before, is 
the noiseless Berry phase and whose variance is 



J 



^2 c)2 



12 



TT COS i9o sin 
TB 



-TT 



(r2+w2)2 



rr 



r2 



2(7l 



TT sin 
TB 



TT -I 



-TT-\ 



r2 



(15) 



This expression has an interesting limiting value when 
Ti <^ oj and lo. To first order in TT: 



2 /7rcos??osin7?o 
^""A B- 



-T12T 



-2at 



TT sin t?( 
B 
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(16) 

which 



We see that the leading term is o-| ^il^i^ 

tends to zero for little i9o- 

When r the fluctuating field has time enough 

to make many uncorrelated oscillations during the cyclic 
evolution. In this case the effect of the fluctuations aver- 
ages out and in the limit (TT)^^ ^ do not give contri- 
bution to the variance which tends to zero: 



= 2(7 



12 



vr cos t?o sin 'do 
B 



vr sin'^ z9c 
B 



1 



T3T 
(17) 



This is to be compared to the dynamical phase which 
grows linearly in T. This different behavior is due to 
the fact that while Berry phase corrections are propor- 
tional to 1/T J Kdt, corrections to the dynamical phase 
are proportional to / Kdt. For a OU process the variance 
of the integral grows linearly for times long compared to 
the autocorrelation time of the field. This is analogous 
to the variance of the position of a brownian particle. 

Until now we concentrated only in the geometrical 
phase. However during an adiabatic cyclic evolution the 
eigenstates acquire both the dynamical and geometrical 
phase. It is known that the dynamical phase d is propor- 
tional to the modulus of the magnetic field. This means 
that the dynamical phase becomes a stochastic processes 
like the Berry phase. We can write S in terms of the fields 
B and K as we did for the geometrical phase: 



BK 



B 



-dt 



(18) 



where — BT Note that expression H18|l is similar to 
(|14|l for Berry phase. The difference is that while 7b 
comes from an integral in parameter space, 5 from an 
integral in the time domain. For instance this means that 
if we double time T, 73 scales with T~^ while the domain 
of integration of 5 doubles. As we will see this is crucial 
for the different role of the two phases in dephasing. 

Following the same steps as for the Berry phase it is 
possible to demonstrate that (5 is a stochastic processes 



with a gaussian distribution. Now we analyze the effect 
of noise on the coherence of a system, in other words de- 
phasing. Suppose we prepare the system in a state which 
is a superposition of the two eigenstate of the Hamilto- 



lV'>=a|T)+foU) 



(19) 



After a slowly cyclic evolution the eigenstates have ac- 
quired both the dynamical and geometrical phases and 
the final state is: 



(20) 



where a = + 5 is the total phase. In the presence 
of noise this phase is a random variable with a gaussian 
distribution P{a) then actually the system at the end of 
the evolution is in a mixed state. 

It can be described by the density operator which is 
given by the expression: 



p= / W){i>'\P{a)da 



(21) 



We want to stress that P{a) ^ P{'^b)P{S), i.e. dynami- 
cal and geometrical phases are not independent processes 
because both depend on K. 

If we insert ea. (|20() in the definition of p we find that 
the population are unchanged while the coherence are 
shrunk by a factor exp(— 2(7^). In terms of the Bloch 
vector, this means that the z component is unchanged 
while the component parallel to the xy is reduced. This 
is what is called dephasing because the relative phase in 
a superposition is undefined. 

In order to do not lose in generality and to compare the 
dynamical and geometrical phase we have studied the two 
case together. We have found the probability distribution 
P(a) for a which has mean value (a) = (75) + (<5) and 
variance: 



<ji = 2 



'12 



TT cos "iJo sin I?, 



^0 



52 



T 

mi2~^') 



B sin do 
T12T 



(22) 
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TT sin i?o 
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In eq. H22() cr^ is the sum of two terms one coming from 
fluctuation in z direction and one in the xy plane. Each of 
these terms contains a factor in round brackets in which 
we recognize a geometrical term proportional to 1/T, as 
we found in (|15|l and a dynamical one proportional to the 
Bohr frequency. Now because of the adiabaticity condi- 
tion we have that the first is much less than the second. 
As a consequence the main contribution to dephasing 
has dynamical rather than geometrical origin. This does 
do not mean of course that in a system in which Berry 
phase emerges dephasing is less than in a system in which 
it does not. What we have demonstrated is that fluctu- 
ations in Berry phase do not contribute considerably to 
dephasing. 

Another relevant aspect to stress is that our calcula- 
tion was performed under the assumptions of a constant 
circular precession, however our results are independent 
of the specific path executed by the magnetic field as long 
as the adiabatic approximation is valid. 

It is worth noting that dephasing is not the only deco- 
herence source in our system. The Bloch vector in fact 
does not return to its initial position since the magnetic 
field does not. To calculate the correct Bloch vector we 
have to average the final positions. However the contri- 
butions from this effect are proportional to {Ki/B)'^ and 
so we can neglect this effect at our level of approximation. 

In this paper we have calculated the distribution of 
Berry phase in the presence of classical noise. Assum- 
ing a OU process for noise we found that under the as- 
sumption of small fluctuation Berry phase is a gaussian 



variable. We have calculated its mean value and the vari- 
ance and we found that the variance diminishes as 1/T. 
This is to be compared to the variance of the dynami- 
cal phase which grows linearly with T. This shows that 
adiabaticity and the geometrical aspects of Berry phase 
reduce fluctuations. This is what was expected but never 
demonstrated. Another aspect that makes Berry phase 
more robust than dynamical phase is that of dephasing. 
We showed that geometrical dephasing is much less than 
the dynamical one. This is due mainly to adiabaticity. 
This means that probably quantum gates based upon 
geometrical phase are more resistant. We would like to 
point out that after we submitted our manuscript a pa- 
per has appeared |0 in which the noise is treated fully 
quantum mechanically and similar conclusions have been 
drawn (but in a completely different setting) about the 
invariant of the geometric phase. These two results to- 
gether complement each other. 
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